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;h ' Using the ZEUS detector at HERA, we have studied the reaction e^p — > e'^X for 

Q^ > 5000 GeV^ with a 20.1 pb~^ data sample collected during the years 1994 to 
1996. For Q^ below 15000 GeV^, the data are in good agreement with Standard 
Model expectations. For Q^ > 35000 GeV^, two events are observed while 0.145 it 
0.013 events are expected. A statistical analysis of a large ensemble of simulated 
Standard Model experiments indicates that with probability 6.0%, an excess at 
least as unlikely as that observed would occur above some Q^ cut. For x > 0.55 
and y > 0.25, four events are observed where 0.91 it 0.08 events are expected. 
A statistical analysis of the two-dimensional distribution of the events in x and y 
yields a probability of 0.72% for the region x > 0.55 and y > 0.25 and a probability 
of 7.8% for the entire Q^ > 5000 GeV^ data sample. The observed excess above 
Standard Model expectations is particularly interesting because it occurs in a 
previously unexplored kinematic region. 
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1 Introduction 

Deep-inelastic scattering (DIS) of leptons on nucleons has been an important tool for 
understanding nucleon structure and many elements of the Standard Model, including 
both the electroweak interaction and quantum chromodynamics (QCD). At the HERA 
collider, DIS processes are being studied at a center of mass energy y/s = 300 GeV and 
at Q^ (the negative of the square of the four-momentum transfer) exceeding the squares 
of the weak vector boson masses. In this regime, lepton-nucleon scattering allows unique 
and sensitive tests of the Standard Model as well as of certain extensions to it |I[]. 

This paper presents results from e~^p running with the ZEUS detector during the years 
1994 to 1996, at proton and positron beam energies of Ep = 820 GeV and Ee = 27.5 GeV. 
With the integrated luminosity of 20.1 pb~^ collected in this period, it has become possible 
to study the reaction e~^p -^ e^X in the region where the expected DIS cross section is 
in the subpicobarn range. This region of high Q^ and x (the Bjorken scaling variable) 
has never before been explored. The above reaction is understood to be a positron- 
quark collision with center-of-mass energy y/xs. Initial cross section measurements by 
the ZEUS and HI |Q collaborations are in good agreement with Standard Model 
expectations for Q^ up to about 10"^GeV^. In this paper, we report on a more sensitive 
search for deviations from Standard Model predictions in the region Q^ > 5000 GeV^. 



2 Neutral Current Deep— Inelastic Scattering 

The reaction studied is: 

e+ + p^e+ + X (1) 

where X represents the final state hadronic system. In the high Q^ regime, the Standard 
Model neutral current (NC) cross section for (0) depends on well-measured electroweak 
parameters and on the parton densities in the proton. Though the latter have not yet been 
measured at high Q'^, perturbative Quantum Chromodynamics (pQCD) predicts their 
values through evolution from high-precision measurements made at lower Q^ values. 

The Born cross section |^ for the NC DIS reaction (|l|) with unpolarized beams isQ 

'^''' ^™'{n(y)^2(a:,Q')-r-(y)^-^3(a:,Q')} , (2) 



dx dQ^ xQ^ 

where a is the electromagnetic coupling. The cross section is given in terms of Q^ and the 
DIS scaling variables x and y = Q'^/sx. In the region of large x and Q^ studied here, the 
parity-violating xJ-'^ term substantially reduces the e~^p cross section, while increasing the 
cross section for e~p scattering (where the second term has positive sign). The explicit 
^/-dependence, which is due to the helicity dependence of electroweak interactions, is 
contained in the functions 
Y±{y) = l±{l-y)\ (3) 

^We neglect the contribution to the cross section (0) of the longitudinal structure function, JF^, which 
we estimate from pQCD and the parton densitiespl to be less than 1% in the kinematic range under 
study. 



while the dependence on the quark structure of the proton, and on the Z'^ propagator is 
absorbed in the (positive) structure functions: 

( j-2ix,Q') \_^ ^ ( cim[qix,Q') + qix,Q')]\ 

written in terms of the quark densities in the proton (g = u, d, c, s, t, b) and the 
corresponding antiquark densities q. For e~^p scattering, the Q^-dependent coefficient 
functions, Cf and Cf, are given by: 

CI{Q')= el -2e,v,v^xz + {vl + al){vl + al)xl .^. 

CKQ'^) = -2egayaeXz + i'2vgag){2veae)xl 

with 

In eqs. |^ and ^ Mz is the Z^ mass, e^ is the quark charge in units of the positron charge, 
Vg = {T^q — 2eq sin^ 6^,) and a^ = T^g are the vector and axial vector couplings of the quark 
to the Z°, Ve and Oe are the corresponding electron couplings, 9^, is the weak mixing angle, 
and T3 is the third component of the weak isospin. All relevant electroweak parameters 
have been measured to high precision 0. 

The QCD-evolved structure functions of equation (^), evaluated at a given x at high 
Q^, depend on quark and gluon densities in the proton measured at lower values of Q^ 
and higher values of x. At high x, u quarks give the dominant contribution to the cross 
section because they have the largest density [Q and because e^ = 2/3. In addition, the 
antiquark (g) density is small 0. 

Uncertainties in the Born-level e^p DIS cross section predictions in this region of high x 
and Q^ are estimated to be about 6.5% (see Section ^, mainly due to uncertainties in 
the evolved quark densities. 

It should be noted that an anomalously high cross section for the production of jets with 
high transverse energy in pp collisions, as recently reported by the CDF collaboration 



T0| , can be explained by adjusting the gluon density in the proton [|TT] (which raises 
the rate of gluon-quark collisions at high x), rather than by adjusting quark densities. 
This variation of the gluon density, however, has only a small effect on the cross section 
predictions relevant to this paper (see Section H). 

3 ZEUS Detector and Monte— Carlo Simulation 
3.1 Experimental Setup 



A description of the ZEUS detector can be found in references ||1^ , [13[. The primary 
components used in this analysis were the compensating uranium-scintillator calorimeter, 
the central tracking detector, and the luminosity detector. 



The calorimeter \Ti\ is divided into three parts, forward (FCAL) covering the polar angle^ 



interval 2.6°<^<37°, barrel (BCAL: 37°<^<129°) and rear (RCAL: 129°<^<176.1°). The 
calorimeters are subdivided into towers which each subtend solid angles from 0.006 to 
0.04 steradians. Each tower is longitudinally segmented into an electromagnetic (EMC) 
section and two hadronic (HAC) sections (one in RCAL). Each HAC section consists of 
a single cell, while the EMC section of each tower is further subdivided transversely into 
four cells (two in RCAL). In test beam conditions, for particle energies up to 120 GeV, 



energy resolutions of aE/E=18%/JE{GeV) for electrons and aE/E=35%/JE{GeV) for 
hadrons have been measured. The cell-to-cell variations in the energy calibration are 
approximately 2% for the EMC cells and 3% for HAC cells. The FCAL and BCAL 
energy scales are presently understood to an accuracy of 3%. The time resolution is 
below 1 ns for energy deposits greater than 4.5 GeV. The impact point of the scattered 
positron at the calorimeter, determined using pulse height sharing, has a resolution of 
about 1 cm. 

In the physics analysis, only those calorimeter cells with energy deposits above thresholds 
of 60 MeV and 1 10 MeV for EMC and HAC cells respectively were used. 

The central tracking chamber (CTD) [|T3] operates in a 1.43 T solenoidal magnetic field. 
It is a drift chamber consisting of 72 cylindrical layers, organized into 9 superlayers. A 
momentum measurement requires a track to pass through at least two superlayers, cor- 
responding to a polar angle region of 15°<^<164°. The transverse momentum resolution 
is (T{pt)/pt = [0.005pi(GeV)] © 0.016 for full length tracks. For full length tracks with 
momenta p > 5 GeV the vertex resolution is 0.1cm in the transverse plane and 0.4 cm 
along Z. 



Events were filtered online by a three-level trigger system |T^ . The trigger criteria used in 
this analysis relied primarily on the energies measured in the calorimeter. The first level 
trigger decision was based on electromagnetic energy and total transverse energy {Ef). 
The second level trigger rejected backgrounds (mostly p-gas interactions) for which the 
calorimeter timing was inconsistent with an ep interaction. In addition, the second level 
trigger applied increased Ef thresholds and also required a minimum value oi E — pz (see 
Section ^, where E and pz are the summed energy and Z-component of the momentum 
measured in the calorimeter. The third level trigger applied more stringent timing cuts 
as well as increased energy and E — pz thresholds. In all cases, the requirements were 
less stringent than those imposed by the offline event selection. 

The luminosity was measured by the rate of high energy photons from the process 



ep —>■ ep7 detected in a lead-scintillator calorimeter [IG] located a.t Z = —107m. The 



uncertainty associated with luminosity measurements is addressed in section |. 

3.2 Monte Carlo Simulation 

NC DIS events were simulated using the meps option of lepto |T^ interfaced to HER- 
ACLES [|l|] via DJANGO [|l9l and the MRSA parton distribution set [|3l- The event sim- 



ulation included electroweak radiative corrections, leading order QCD effects and parton 



showers. Hadronization was simulated with JETSETf21 



^Thc right-handed ZEUS coordinate system is centered on the nominal interaction point (Z = 0) 
and defined with the Z axis pointing in the proton beam direction, and the horizontal X axis pointing 
towards the center of HERA. 



Large samples of simulated photoproduction events [^] were used for background studies. 
Samples of both direct and resolved photoproduction events (including the production of 
cc and 66 pairs) were generated using both pythia |^ and herwig |]23|. Direct and 



resolved photoproduction of events with prompt photons were simulated with herwig. 



Production of W and Z bosons was studied using the epvec [^] generator. Finally, the 



processes 77 — > e+e and 77 -^ t^t were simulated using zlpair |25| 



All MC events were passed through a geant [|26[ based simulation of the ZEUS detector 
and trigger, and analyzed with the same reconstruction and offline selection procedures 
as the data. 



4 Positron Identification and Event Kinematics 

A key signature of high Q^ e^p — > e^X events is an isolated high transverse momentum 
positron. In order to identify and reconstruct this positron, while rejecting events in 
which other final state particles mimic a positron, an algorithm was used which combines 
calorimeter and CTD information. 

In a first step, the calorimeter cells are clustered by joining each cell to the highest energy 
cell among its adjacent neighbours. All clusters are evaluated as positron candidates. The 
cluster energy, -Eciu, is the sum of the cell energies belonging to the cluster. The cluster 
angle, 6'ciu, is set equal to the polar angle obtained by joining the energy-weighted mean 
position of the cluster with the event vertex obtained from the tracks measured with 
the CTD. For candidates with polar angleQ within the CTD acceptance (^ciu > 17.2°), 
a matching track is required. A track is considered to match if the distance of closest 
approach (DCA) between the extrapolation of the track into the calorimeter and the 
position of the cluster center is less than 10 cm, where the r. m. s. resolution on the DCA 
is 1.8 cm. 

In the second step, several quantities, ^j, are calculated for each positron candidate: the 
fraction of the cluster energy in the HAC sections of the calorimeter, the parameters 
related to lateral energy profiles, and the total energy (-Econe) in all calorimeter cells 
not associated with the cluster but lying within an 77, cj) (pseudorapidity,azimuth) cone of 
radius 0.8 centered on the cluster. If a matching track is present, we also evaluate the 
polar and azimuthal angle differences between the track and the cluster position, and the 
quantity 1/-Eciu ~ l/-ftrk) where Ptrk is the track momentum. 

Finally, we transform each ^j into a quality factor Q{ii). Candidates are accepted as 
positrons if the product of the Q{^i) exceeds a threshold determined from Monte Carlo 
studies. The efficiency for finding positrons in a neutral current DIS sample with Q^ > 
5000 GeV^ is 91%. In accepted events, the positron energy, E'^, is set equal to the cluster 
energy, Edu, and the positron angle, 9e, is set equal to 6'ciu- The resolution in 9e is typically 
better than 0.3°. 

For each event with an accepted positron, the following global event quantities were 



^We do not consider candidates with 6*0111 > 164° (which are also beyond the CTD acceptance hmit), 
since they correspond to Q^ values below the range of this analysis. 
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calculated from the energy deposits in the calorimeter: 



Pt 

E-pz 


/ \2 / \2 

= \ Y.fx + Y.Py 

i 


Et 


i 


(Pt)had 
{E - pz)had 


= \ T.PX + y,fy 



(7) 



where the sums run over all calorimeter cells with energy deposits above threshold and 
the p* are the momenta assigned to each calorimeter cell (calculated assuming zero mass 
with the direction obtained from the cell center and the measured vertex position). The 
primed sums exclude the cells associated with the positron. 

To describe the hadronic system, we use the angle, 7raw5 and energy. Eg, defined as 



cos 7i., 



(Pf)had -{E- pz)Ld 

{Pt)Ld + {E- Pz)had 



and 



E„ 



{P 



t had 



Sm 7i.aw 



Resolution effects and systematic shifts of 7i.aw have been studied with MC simulations. 
The reconstructed 7raw is systematically higher than the generated value by about 2.7°. 
To remove this bias, we compute a corrected value, 7, which depends on 7raw and Oe- The 
r. m. s. resolution of 7 is about 2.5° for x > 0.55 and Q^ > 5000 GeV^. 

In the quark-parton model, for a perfect detector, 7 and Eg are interpreted as the scat- 
tering angle and energy of the massless quark q in the reaction eg — *> eg. 

At a given value of s, the kinematic variables (x, y, and Q"^) can be reconstructed from 
any two of the four measured quantities: E'^, 6e, Eg, and 7. Different combinations have 
been used by the HERA experiments. At high x and Q^ where the calorimeter energy 
resolution functions are narrow, the dominant uncertainties in energy measurements are 
due to systematic effects such as energy loss in inactive material in front of the calorime- 
ter, nonuniformities and nonlinearities in the calorimeter response, longitudinal energy 
leakages, and energy carried away by neutrinos and muons. For the hadronic system, the 
raw measured energies are typically 15% less than the true energies. For positrons, the 
raw measured energies are typically 4% less than the true values. 

We choose the double-angle method |^ because it is least sensitive to uncertainties in 
the energy measurement. In this scheme, the kinematic variables are obtained from O^ 
and 7 as follows: 



K 



2/da 

Q 



sm7 



sin^p 



2 

DA 



Ep (1 — C0S7) (1 — cos6'e) 

sin 6^6 (1 — cos 7) 
sin 7 + sin 9e — sin(7 + 9e) 

S 3;da 2/da • 



(9) 
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For y > 0.25 and x > 0.45, the resolution in x^^ is 9%; it improves to 6% for y > 0.5. 
The resolution in Q^^ is typically 5% at large x and y. 

For selected events with high x and high Q^ we also present the kinematic variables 
calculated from the scattered positron energy E'^ and angle 9e using the equations: 



Ee Elil + cosOe 






Ep 2Ee- E'^il-cosOe 



ye = l-^(l-COS^e), (10) 



2E, 



e 



C^g S Xg ye • 

We apply a test-beam based correction to E'^ to account for energy loss in inactive material 
and nonuniformities of the calorimeter response. 



5 Event Selection 

Important characteristics of reaction (|l|) that distinguish it from background processes 
include (i) the presence of an energetic isolated positron, (ii) pt balance, and (iii) E—pz ~ 
2Ee = 55 GeV. In addition, at large Q^, the transverse energy Et typically exceeds 
100 GeV. 



About 10^ events were accepted by the trigger requirements described in section |3.1| . The 
offline event selection criteria are described below. 

• E -pz 

The net E — pz as measured in the calorimeter is required to be in the range 
40 GeV < E-pz < 70 GeV (44 GeV < E - pz < 70 GeV) for 6^ > 17.2° {9^ < 
17.2°). The lower cut rejects backgrounds such as photoproduction or e~^p -^ e+X 
events with a hard initial state photon, for which energy escapes through the rear 
beam hole (see below). The 70 GeV cut removes a small number of events with a 
misreconstructed vertex position. 

• Longitudinal vertex position 

The event vertex reconstructed from CTD tracks must have a Z position {Z^tx) 
within 50 cm of the nominal interaction point. The Z^tx distribution of the data is 
roughly Gaussian with {Z^tx) = —2 cm. The r. m. s. spread in Z^tx, 12 cm, is largely 
due to the length of the proton beam bunches. 

• Positron requirements 

An isolated positron candidate with energy E'^ > 20 GeV and -£"00110 < 5 GeV must 
be found by the algorithm described in section ^ Additional requirements depend 
on the polar angle of the positron: 

For 6e > 17.2°, where the positron candidates are within the CTD acceptance, 
a matching track with momentum above 2 GeV is required. 
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For Of. < 17.2°, where the positron either misses the CTD altogether or is on 
the edge of the CTD acceptance, the number of fake positron candidates is 
large. These have a sharply falling transverse momentum spectrum. To reduce 
this background, we require positron candidates in this angular range to have 
transverse momenta above 30 GeV. 

To remove Compton scattering events [ep — > e'jX), we reject any event which has 
two isolated electromagnetic clusters in the calorimeter, each with Edu > 8 GeV and 
£;cone < 2 GeV. 

• Momentum transfer 

We require Ql^ to exceed 5000 GeV^. 

The overall selection efficiency, estimated using Monte Carlo NC events generated with 
Q^ > 5000 GeV^, is 81%. For the 191 events which pass all cuts, the mean measured 
E — pz is 51.9 GeV with an r. m. s. width of 4.2 GeV, in good agreement with the Monte 
Carlo e^p NC simulation which predicts a mean of 51.8 GeV and an r. m. s. of 4.0 GeV. 
While no cut was applied to the net transverse momentum (pt), the surviving events have 
a mean pt of 7.5 GeV, again in good agreement with the e~^p NC Monte Carlo prediction 
of 7.1 GeV. 



6 Data and Expectations at Large x and Q^ 

Figure |l] shows the distribution in the (xda,1/da) plane of the 191 events satisfying the 
selection criteria. In Table [^, the numbers of observed events are compared with the 
Standard Model expectations in bins of Xda and ?/da- In general, the agreement between 
the data and the Standard Model expectations is good. However, five events, in four 
(3;da,?/da) bins occur at high x^a and Q^^ where the expected numbers of events are 
small. Four lie in the region x^a > 0.55 and ?/da > 0.25, while the fifth has Xda = 0.48 
and a very high Q"^^. These five events are selected for more detailed discussion below. 

Figures and ^ show the x^a (for y^j^ > 0.25) and Q^a distributions of the final event 
sample. In both figures, the e~^p NC prediction for the same integrated luminosity is 
superimposed as a solid histogram. Again, the agreement with the Standard Model is 
good at lower values of x-qa and Q"^/^, but an excess is observed at high Xda and at high 

Table |^ shows the kinematic variables, before applying the corrections discussed in section 
H, associated with the five selected events. Included are the uncorrected values of x^a, 
j/da, and Q^A (calculated using 7raw) as well as the corrected value of 7. Table ^ gives 
the kinematic variables and their estimated uncertainties obtained using the double-angle 
and electron methods. The uncertainties have been estimated from the resolutions in 7 
and 6e, as well as estimates of the systematic uncertainty in the 7-correction procedure 
discussed in Section ^. The quoted r. m. s. errors on the electron variables include the 
uncertainty in 6^, the calorimeter energy resolution, the uncertainty associated with the 
calorimeter nonlinearity, and the uncertainty on corrections applied for inactive material 
and nonuniformities. Though 6^ is used in both the DA and electron methods, it makes 
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only a small contribution to each error. Hence the errors on the two measurements are 
essentially independent. 

All events listed in Tables ^ and|^, except the first, have a track matching the electromag- 
netic shower of the scattered positron in the calorimeter. In these events, the positron 
track momentum is consistent with the calorimeter energy within measurement error^. 
The first event (ll-Oct-94) has a positron candidate at too small an angle to produce an 
observable track in the CTD.Q We show event displays of the first two events in Figs. ^ 
and |. 

The five events have clean, well-identified and isolated positrons and jets in the final state. 
None lie close to any of the selection cuts described in the previous section. For these 
events, the scattering angles and energies of the final state positrons and jets are measured 
with good precision, making it unlikely that resolution smearing has moved any of these 
events from low Q^ to the measured Q^j^. 

Initial state radiation (ISR) from the incoming positron, where the radiated photon es- 
capes through the rear beam hole is a possible source of uncertainty in the determination 
of the event kinematics. Since ISR affects the DA and electron variables differently, it is 
possible to estimate the energy E^ of the radiated photon. For each of the five events, 
E^ is consistent with zero within resolution and the measured values oi E — pz limit 
E^<3 GeV. 



7 Background Estimation 

Potential backgrounds to e'^p DIS events at large x and y are those processes which yield 
an isolated positron or electron of high transverse energy, or a photon or 7r° which could 
be misidentified as a scattered positron. The latter event class contributes predominantly 
to the background of events in which the positron is very forward {6e < 17.2°) and no 

track information is available for the positron candidate (e.g. the first event in Tables @ 
and|^). At larger angles, photon conversions in inactive material between the interaction 
point and the CTD can also mimic positron candidates with matching tracks, but this 
effect, which is included in the detector simulation, is much smaller. 

In the following, we describe the physical processes studied as possible sources of back- 
ground. Limits are quoted at 90% confidence level. 

• Prompt photon photoproduction {'-fp — > 7X) has been studied using HERWIG. We 
generated an event sample with the final state photon transverse momentum ex- 
ceeding 20 GeV. The cross section is 1.6 pb, of which 86% (14%) is due to direct 
(resolved) photoproduction. The observed cross section due to this process in the 
region x^a > 0.45 and y^A > 0.25 is 0.28 fb (0.006 events). 

• Photoproduction of high Et jets can contribute to the background if a jet is misiden- 
tified as a positron. Using herwig, we have generated event samples for both direct 



* It should be noted that the positron energies in table g are so large that the tracking error does not 
allow an unambiguous determination of the particle charge. 

^There are hits in the innermost layer of the CTD, aligned in azimuth with this positron candidate. 
However, the hits are too few to qualify as a track according to our standard criteria. 
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and resolved processes which include heavy quark production and decay. In these 
samples, no event satisfies the selection criteria for Xda > 0.45 and ^da > 0.25, 
providing an upper limit of 1.8 fb (0.04 events). 

• QED Compton scattering [ep —>■ e'jX) could produce background if one of the 
electromagnetic showers is not recognized as such. Monte Carlo studies show that 
this probability is negligible, with an upper limit on the contribution to the observed 
cross section of 0.2 fb (0.004 events). 

• Two photon production of lepton pairs (77 — »• ££) was studied using ZLPAIR. No 
events from the process 77 — > e~^e~ were found after the selections. For 77 — »• t~^t~ , 
where one r decays via r -—>■ eu, the quantity E—pz as well as the electron transverse 
energy are typically much lower than for high Q^ NC events. We obtain the upper 
limit on the contribution to the observed cross section of 0.1 fb (0.002 events). 

• Leptonic decays of W bosons have been studied using a Monte Carlo sample gen- 
erated with EPVEC. The total cross section for production of W^ bosons and their 
subsequent decay via W -^ tv^ is approximately 0.1 pb. The final state contains a 
(anti)neutrino with high transverse momentum (of order 40GeV), which typically 
results in large missing E — pz (as well as pt). We estimate the accepted cross 
section for this process to be less than 0.5 fb (0.01 events). Decays of the neutral 
boson, Z^ — > e~^e~ , are rejected by the cut on two electromagnetic clusters and are 
expected to contribute a negligible background. 

The estimated cross sections from these background sources are listed in Table ^ along 
with the e~^p NC cross section. The backgrounds are much smaller than the DIS signal 
in the region of interest, and are neglected. 



8 Uncertainties of the Standard Model Predictions 

The predicted numbers of e~^p NC DIS events depend on (i) the measured luminosity, (ii) 
the electroweak parameters, (iii) electroweak radiative corrections, mainly due to initial 
state radiation (ISR), (iv) the quark densities in the relevant region of x and Q^ and (v) 
the Monte Carlo simulation of the detector. We now discuss the precision to which these 
quantities are known and describe the studies performed to determine the uncertainties 
of the predictions. 

• Luminosity measurement 

The luminosity is measured to a precision of about 1.5 % using the ZEUS luminosity 
monitor. The recent 1996 running period has a larger uncertainty due to effects from 
beam satellite bunches. Also, the offline calibration of the luminosity detector is 
not yet finalized. Including these uncertainties from recent data, the uncertainty for 
the full data sample is 2.3%. 

• Electroweak parameters 

The relevant electroweak parameters have been measured to high accuracy and 
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contribute a small uncertainty in the predicted cross section over the HERA kine- 



matic range [^. The HERACLES program calculates NC DIS cross sections to first 
order using input values for the Fermi constant G^, M^, the top mass rrit, and the 
Higgsmass. Varying M^ = 91. 187 ± 0.007 GeV and mi = 180 ± 12 GeV within their 
experimental errors [^ changes the predicted cross section in the kinematic range 
reported in this paper by only 0.25%. 

Radiative corrections 

The program HECTOR p9|, which includes the effects of second order QED radiative 
corrections was used to check the cross sections computed using HERACLES. The 
differences were found to be about 1.5% for the integrated cross sections in the 
region x^a > 0.5 and ycA > 0.25. 

The luminosity monitor records data for all triggered events, and so measures di- 
rectly, with an acceptance of about 30%, the ISR spectrum for accepted events. The 
experimental data are in quantitative agreement with the ISR spectrum calculated 
for the accepted sample. 

Corrections due to initial state radiation convoluted with the experimental resolu- 
tion, based on studies |^ made for lower values of x, produce uncertainties of less 



than 2% in the accepted cross sections. This number is used as the estimate of the 
uncertainty due to radiative corrections. 

Structure functions 

The least well known inputs to the predicted cross section in equation ^ are the 
structure functions. To estimate the uncertainty associated with parton densities, we 
performed a NLO QCD fit to fixed-target F2 lepton-proton data (with x > 0.1) from 
the NMC m, SLAG ||], and BGDMS |3| collaborations and xF^ and q/xFs resuhs 
from the GGFR collaboration |^. A complete treatment of statistical and correlated 
experimental systematic errors was included in the fit. The results of the fit are 
consistent with the MRSA |^ and GTEQ3 |^ parton density parameterizations 



up to Q2 of 5 X lO^GeV^. 

The fit was used to estimate the two largest uncertainties due to the structure 
functions: the experimental uncertainties and the uncertainty of the quark-gluon 
coupling, as, used in the evolution to higher Q^. The effects of experimental un- 
certainties in the fixed-target data result in a ±6.2% uncertainty in the integrated 
cross section at HERA for x > 0.5 and y > 0.25. The uncertainty due to a^ was 
estimated by varying the value of as{Mz) used in the QGD evolution from 0.113 to 
0.123, which produces an uncertainty of ±1.9%. From the above studies, we take 
the overall uncertainty in the cross section due to structure function uncertainties 
to be ±6.5% over the kinematic range of interest. 

Other sources of uncertainty in the structure functions were found to be small. 
Changing the strange quark fraction in the QCD fit from 10% to 30% produced less 
than 0.1% change in the predicted cross section. Removing BGDMS data from the 
fit produced a change of only 1.7%. Removing data with W"^ = sy{l — x) between 
10 and 25 GeV^ had no significant effect. Since the contribution of charm to the 
cross section for x > 0.5 and y > 0.25 is 0.5%, uncertainties in the charm quark 
mass and the charm evolution renormalization scale can be safely neglected. 
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As a cross check, the uncertainty of 6.5% was compared to the differences in cross 
section predicted by various parton density parameterizations. For example, a com- 
parison of integrated cross sections predicted by the MRS A, CTEQ3, and GRV94 



35[| parameterizations produces an r. m. s. of 2%. A comparison of the CTEQ4 HJ 
parameterization |11| (which was tuned to the CDF high Et jet cross section |p!0| ) 
with the nominal CTEQ4 parameterization produced an increase in cross section 
of only 1.9%, demonstrating the small effect at HERA of a larger gluon density 
at high X. Finally, a crude estimate of the contributions from QCD corrections at 
higher than NLO can be estimated by comparing the cross sections predicted by the 
GRV94 LO and NLO parameterizations, which produced a cross section difference 
of only 1%. 

Table |^ summarizes the structure function uncertainties as well as the cross checks 
which were performed. 

• Detector simulation 

To estimate the uncertainties in the expected event yields due to possible inaccu- 
racies in the detector simulation, we made several modifications to the simulation 
to reflect uncertainties in the overall calorimeter energy scale and in the simulation 
of the calorimeter and CTD response to positrons. The FCAL and BCAL energy 
scales were separately varied by ±3%, our present estimate of this uncertainty. Each 
of the seven measured quantities used in the positron identification algorithm was 
varied by an amount consistent with the differences between the data and the nom- 
inal simulation. For the region Xqa > 0.55 and |/da > 0.25, the resulting uncertainty 
in the expected number of events is 4.4%. 

We conclude that at the large x and Q^ values discussed in this paper the overall uncer- 
tainty of the number of events predicted within the Standard Model is 8.4%. 



9 Comparison of Data with Standard Model and Sig- 
nificance of Excess 

Table [I| compares the data with the e~^p -^ e^X expectations in bins of Xda and |/da 
for Q^DA > 5000 GeV^. There is very good agreement over the entire plane, except in 
the region of high x^a and ?/da- The numbers of observed and expected events above 
various Qda thresholds are given in table ||. The data agree well with the Standard Model 
predictions up to Qda of 1-5 x lO'^GeV^. 

Fig. ^ shows the number of events with Q"^^ > Qf^j^ as a function of Qda- Figure ^ 
shows the number of events with |/da > 0.25 and with Xda > a^DA^ ^s a function of x*a- 
On each of the two plots, the e~^p NC DIS Monte Carlo expectation is shown as a dotted 
line. 

We define the Poisson probability corresponding to the event numbers in Fig. ^a as 

oo n 

'PiQll)= E ^e~' (11) 
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where A^obs is the number of observed events with Q^a > Qda) ^^^ A^ i^ ^^^ number 
of events expected from NC DIS in the same region. In Fig. §3 V{Q^jJ is shown as a 
function of Q^a- The minimum probabihty of V{Q'^j^ = 0.39% (corresponding to 2.7 
Gaussian standard deviations) occurs at Q^^ = 3.75 x 10^ GeV^ where two events are 
observed while 0.091 ± 0.010 are expected. If the expected number of events is increased 
by its error, V{Q^j^) increases to 0.47%. 

We have performed a similar analysis of the Xda spectrum in the region ^oa > 0.25. The 
probability V{x^jJ is shown as a function of x*^ in Fig. |^. Here the minimum value 
'P(x*^) = 0.60% (corresponding to 2.5 Gaussian standard deviations) occurs at x*^ = 0.57 
where four events are observed and 0.71 ± 0.06 are expected. If the expected number of 
events is increased by its error, 'P(x*^) increases to 0.79%. The corresponding results for 
different ?/da cuts appear in Table |^. 

To gauge the significance of these probabilities, one must consider that it is possible 
to observe a statistical fluctuation above any Q^l or x*^ within the region studied. 
We generated a large ensemble of simulated experiments according to Standard Model 
assumptions, each with a luminosity of 20.1 pb~^ and asked how often an experiment 
would have a probability 'P(QdI) < 0.39% for any Q'^j^. The resulting probability to 
find such a fluctuation was 6.0%. Similarly, we determined that the probability for an 
experiment to have V (x*^) < 0.60% in the region ^da > 0.25 for any Xda was 7.2%. The 
same analysis was applied for other y^^ cuts and the results appear in Table 0. 

Finally, we have performed a statistical analysis which computes a probability for the two- 
dimensional distribution of the events in the (xda,2/da) plane (with Q^^ > 5000 GeV^). 
Here the data from each simulated experiment were binned as in Table |1|. Over a given 
region TZ of the (x^a, 2/da) plane, which is defined as a subset of the bins shown in Table |l|, 
we compute the likelihood for a given experiment as 

where iVj is the number of events observed and fii is the number of events expected in bin 
i. For region TZ, we denote by C^^^ the value of C-n obtained from the data. 

Using the ensemble of simulated experiments, we determined the probability that C-n < 
Cj^^ for several choices of the region TZ. If 71 is the entire (xda, Z/da) plane, the probability 
that Cti < C^^ is 7.8%. liTZ consists of the entire (xDA,yDA) plane, except for x^a > 0.55 
and j/da > 0.25, the probability that C-ji < C^^ is 50.2%, indicating that the data are in 
good agreement with the Standard Model in this region. In contrast, the probability that 
Cii < C^' in the region 7^ defined by Xda > 0.55 and |/da > 0.25 is 0.72%. 



10 Conclusions 

Using the ZEUS detector at HERA, we have studied the reaction e~^p -^ e~^X for Q^ > 
5000 GeV^ with a 20.1 pb^^ data sample collected during the years 1994 to 1996. 

For Q^ below 15000 GeV^, the data are in good agreement with Standard Model expec- 
tations. For Q^ > 35000 GeV^, two events are observed while 0.145 ± 0.013 events are 
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expected. A statistical analysis of a large ensemble of simulated Standard Model experi- 
ments indicates that with probability 6.0%, an excess at least as unlikely as that observed 
would occur above some Q^ cut. 

For X > 0.55 and y > 0.25, four events are observed where 0.91 dz 0.08 events are expected. 
A statistical analysis which assigns a probability to the two-dimensional distribution of 
the events in x and y yields a probability of 0.72% for the region x > 0.55 and y > 0.25 
and a probability of 7.8% for the entire Q^ > 5000 GeV^ data sample. 

The observed excess above Standard Model expectations is particularly interesting be- 
cause it occurs in a previously unexplored kinematic region. 
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0.65 


0.65 
0.75 


0.75 
0.85 


0.85 
0.95 


0.95 < t/DA < 1-00 


0.15 


0.015 


0.033 


0.013 


0.0055 


0.0015 


0.0012 






0.85 < t/DA < 0.95 


8.8 
9 


1.2 
3 


0.32 


0.10 


0.028 
1 


0.01 


0.0034 






0.75 < 1/da < 0.85 


12 
16 


2.5 
4 


0.50 
1 


0.15 


0.050 


0.011 


0.0039 






0.65 < t/DA < 0.75 


13 
10 


3.7 
3 


0.86 


0.26 


0.082 


0.022 


0.0054 
1 


0.0020 




0.55 < t/DA < 0.65 


15 
12 


6.1 
3 


1.65 
3 


0.46 
1 


0.15 


0.046 


0.0090 


0.0024 




0.45 < t/DA < 0.55 


12 
6 


11 
13 


2.5 
1 


0.85 


0.28 


0.084 
1 


0.0208 


0.0032 




0.35 < 1/da < 0.45 


4.6 
3 


18 
17 


5.5 
6 


1.75 


0.52 


0.16 


0.0403 


0.0093 




0.25 < t/DA < 0.35 




18 
23 


11 
6 


3.74 
7 


1.19 
1 


0.34 
2 


0.1104 


0.0175 


0.0066 


0.15 < t/DA < 0.25 




2.2 
1 


14 
15 


9.6 
10 


3.32 
3 


1.2 


0.2784 
1 


0.0717 


0.0077 


0.05 < i/da < 0.15 








1.3 
1 


2.14 
3 


1.6 
2 


0.9052 
1 


0.3022 
1 


0.1216 



Table 1: The observed numbers of events in bins of x^^ and y^A (bottom number in 
each box), compared to the expected number of e~^p NC events (top number in each box). 
There are no events observed above x^a = 0.95. 
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ZEUS 1994-1996 



Event Date 


ll-Oct-94 


03-NOV-95 


12-Sep-96 


12-Oct-96 


21-NOV-96 


Et GeV] 


123. 


217. 


193. 


204. 


187. 


Pt GeV] 


8.9 


8.2 


2.9 


2.2 


10.2 


E-pz GeV 


47.8 


53.2 


49.7 


50.2 


49.1 


Eq GeV] 


67.4 


235. 


270. 


151. 


276. 


Traw 


69.0° 


28.1° 


19.9° 


40.7° 


19.7° 


E'^ GeV] 


324. 


220. 


149. 


366. 


134. 


Oe 


11.9° 


27.8° 


39.3° 


15.4° 


41.1° 


(,^DA jraw 

lyDAJraw 
(QL)raw[10^GeV2] 


0.468 
0.868 
3.67 


0.541 
0.503 

2.45 


0.535 

0.330 

1.59 


0.668 
0.733 

4.42 


0.515 
0.316 

1.47 


7 


67.6° 


26.7° 


17.3° 


38.6° 


17.0° 



Table 2: Measured variables for the five events selected as described in the text. The first 
row shows the date the event was acquired. The following rows indicate the quantities 
defined in equations 7 and 8, followed by the energy and angle of the scattered positron. 
The values of x, y, and Q^ calculated from 7i.aw and 6e are shown next. The last row 
shows the 7 angle. 
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ZEUS 1994-1996 



Event Date 


ll-Oct-94 


03-NOV-95 


12-Sep-96 


12-Oct-96 


21-NOV-96 


3^ DA 
"3^DA 


0.480 
0.035 


0.570 
0.029 


0.617 
0.054 


0.709 
0.034 


0.597 
0.053 


1/da 
^I/da 


0.865 
0.008 


0.490 
0.010 


0.299 
0.017 


0.721 
0.008 


0.285 
0.017 


Ql^ [10^ GeV2] 
(5QL lO^GeV^] 


3.75 
0.26 


2.52 
0.07 


1.66 
0.05 


4.61 
0.16 


1.54 
0.04 


(5Xe 


0.525 
0.048 


0.536 
0.048 


0.562 
0.102 


0.605 
0.060 


0.443 
0.063 




0.854 
0.018 


0.505 
0.024 


0.319 
0.039 


0.752 
0.021 


0.350 
0.032 


Ql lO^GeV^] 
bQl [lO^GeV^] 


4.05 
0.31 


2.44 
0.11 


1.62 
0.09 


4.10 
0.30 


1.40 
0.07 



Table 3: Kinematic variables for the five events selected as described in the text. The 
first six lines below the event dates show the double angle values and their estimated 
uncertainties. These include the r. m. s. errors as well as small contributions from the 
uncertainties associated with the correction procedure. The last block of six lines shows 
the kinematic variables reconstructed from the energy and the angle of the positron. 
These latter errors are dominated at present by systematic uncertainties associated with 
the positron energy measurement. 



Background 
Process 


cross sec 
Xda > 0.45 


tion fb 
Xda > 0.55 


7P -^ 7X 
7P — s> dijets 


0.28 
< 1.8 


0.28 
< 1.8 


ep — s> e'^X 


<0.2 


<0.2 


77 -^ ii 


<0.1 


<0.1 


W ^ eu 


<0.5 


<0.5 


Expected NC DIS 


165 


46 



Table 4: Expected cross sections for different background processes in the regions (xda > 
0.45, i/da > 0.25) and (xda > 0.55, |/da > 0.25). The expected numbers of background 
events are obtained by multiplying these cross sections with the integrated luminosity of 
0.02 fb^^. The quoted limits are at 90% CL. Shown for comparison in the last row are 
the cross sections expected for e^p NC events. 
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Systematic errors 


fixed-target experimental uncertainties 


±0.062 


0.113 < a, < 0.123 


±0.019 


overall assumed r. m. s. uncertainty 


±0.065 


Cross checks 


10% < strange fraction < 30% 


< 0.001 


uncertainties in charm evolution 


< 0.005 


GRV94, MRSA, CTEQ3 comparison 


±0.020 


GRV94 NLO versus LO 


±0.010 


High-x gluon (CDF inspired, CTEQ4 HJ) 


±0.019 



Table 5: Relative uncertainties in the integrated cross section for x > 0.5 and y > 
0.25 due to variations in the structure functions. The top two entries represent the two 
dominant contributions to these uncertainties, and so provide the systematic error, shown 
in the third row, which is used in this paper. The remaining entries are cross checks that 
are not independent of the items in the first two rows. 



ZEUS 1994-1996 



Qll GeV2] 


NoUQL > Ql\) 


li 


6fi 


5000 


191 


196.5 


±9.87 


10000 


33 


32.18 


±2.04 


15000 


12 


8.66 


±0.66 


20000 


5 


2.76 


±0.23 


25000 


3 


1.01 


±0.09 


30000 


2 


0.37 


±0.04 


35000 


2 


0.145 


±0.013 



Table 6: The observed and expected numbers of events above various Qda thresholds. 
The first two columns give Q'^j^, the lower limit on Q'^^, and A'^obs, the number of observed 
events with Q^a > Qda- The next two columns give /i, the expected number of events 
with Qda > Qda5 s-^d 6fi, the uncertainty on /x, which includes uncertainties in the cross 
section prediction as well as experimental uncertainties. 
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ZEUS 1994-1996 



yoA range 


' minl^^DAJ 


X* 
•^DA 


A^obs(a;DA > xU 


/i 


-PsM 


I/da > 0.05 


1.61% 


0.708 


4 


0.95 


16.0% 


I/da > 0.15 


2.57% 


0.708 


2 


0.25 


23.0% 


1/da > 0.25 


0.60% 


0.569 


4 


0.71 


7.2% 


I/da > 0.35 


3.38% 


0.708 


1 


0.034 


26.6% 


1/da > 0.45 


1.32% 


0.569 


2 


0.17 


12.7% 


2/da > 0.55 


0.96% 


0.708 


1 


0.010 


9.5% 


I/da > 0.65 


0.50% 


0.708 


1 


0.005 


5.0% 



Table 7: Minimal Poisson probabilities associated with the Xda distributions for different 
I/da cuts. The columns labelled Pminla^DA) and a;*^ give the minimal probability and the 
cut on Xda where it occurs. The next two columns give A'^obs and n, the number of events 
observed and the number expected with Xda > a^oA- The column labelled Psm gives 
the probability that a simulated e~^p Standard Model experiment yields a lower value of 



V„ 



[x: 



than the one observed. All values are for Q^a > 5000 GeV^ 
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0.9 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



X 
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Figure 1: The distribution of the event sample in x-^j^ and ^/da- The hues indicate 
constant values of Ql^ = x^j^yj^j^s for Q^^ = 5000, 10000, 20000 and 40000 GeV^. 
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ZEUS 1994-1996 
y^,>0.25 

qI^>5000 GeV^ 
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Figure 2: The Xda distribution of the observed events with the cuts shown (full dots), 
compared to the Standard Model e~^p NC expectation (histogram). The error bars on the 
data points are obtained from the square root of the number of events in the bin. 
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Figure 3: The Qda distribution of the observed events (full dots), compared to the 
Standard Model e~^p NC expectation (histogram). The error bars on the data points are 
obtained from the square root of the number of events in the bin. 
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ll-Oct-1994 




V 



se^' 







Figure 4: A display of the high Q^ event recorded on ll-Oct-94. The top right part 
shows the ZEUS inner tracking system and the calorimeter. The filled rectangles in the 
calorimeter denote energy deposits which are above the noise thresholds described in 
the text (cf. Section 3.1). The bottom right display shows a projection onto a plane 
perpendicular to the beam axis, where only BCAL energy deposits are shown. The 
left part of the figure shows the calorimeter transverse energy deposits. This display 
demonstrates that the scattered positron is well isolated. 
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03-NOV-1995 



Et (GeV) 




O' -3 



V 



se^' 







Figure 5: A display of the high Q^ event recorded on 03-Nov-95. The description of the 
display is identical to the previous figure. However, for this event the positron polar angle 
de is large enough to be in the CTD acceptance. 
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Figure 6: In (a), the solid line indicates the number of observed events with Q^a > Q'da 
as a function of Q"^^. The dotted line indicates the number of events expected from e^p 
NC DIS with QIj^ > Q'^j^. In (b) is shown the Poisson probability (eqn. |ll|) to observe at 
least as many events as were observed with Q"^^ > Q"^^ as a function of QdI- 
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ZEUS 1994-1996 
y^,>0.25 

qI^>5000 GeV^ 
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Figure 7: In (a), the solid line indicates the number of events observed with y^A > 0.25 
and Xda > a^DA ^^ a function of x* a- The dotted line indicates the number of expected e~^p 
NC DIS events with y^A > 0.25 and Xda > a^DA- ^^ (b) is shown the Poisson probability 
(eqn. |Tl|) to observe at least as many events as were observed with Xqa > a^DA ^^ ^ function 
of X* . 
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